Confined Shocks inside Isolated Liquid Volumes — A New Path of Erosion? 
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The unique confinement of shock waves inside isolated liquid volumes amplifies the density of 
shock-liquid interactions. We investigate this universal principle through an interdisciplinary study 
of shock-induced cavitation inside liquid volumes, isolated in 2 and 3 dimensions. By combining high- 
speed visualizations of ideal water drops realized in microgravity with smoothed particle simulations 
we evidence strong shock-induced cavitation at the focus of the confined shocks. We extend this 
analysis to ground-observations of jets and drops using an analytic model, and argue that cavitation 
caused by trapped shocks offers a distinct mechanism of erosion in high-speed impacts (> 100 m s~^). 
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Shock waves in liquids are a common cause of cavita- 
tion [1-5], in particular when shocks are reflected and 
focussed 0. Famous examples include shock-induced 
cavitation in lithotripsy [tI, [H and the 'white crown' 
on the sea surface following an underwater detonation 
(oMll]. However, little is known about shock-induced 
cavitation inside "isolated" liquid volumes [l2j, which 
are completely bounded by a free surface. The crucial 
feature of such systems is their unique confinement: 
the closed surface acts as a mirror trapping the shock 
and amplifying its local interaction with the fluid. Ex- 
perimental hints for the importance of this mechanism 
for generating cavitation were provided by some of the 
earliest high-speed visualizations of shocked drops [Ref. 
[13I , see reprint in Fig. [1^]. 

In this letter, we study the amplification of shock- 
induced cavitation in liquid volumes isolated in three 
dimensions (3-d), such as drops, and in two dimensions 
(2-d), such as jets. Illustrations of shock-driven cavi- 
tation in both cases are provided in Fig. [H To under- 
stand this cavitation, we perform a systematic experi- 
mental study of shock-induced cavitation inside large, 
spherical water drops. These drops are realized in mi- 
crogravity conditions aboard parabolic flights (Euro- 
pean Space Agency, 42nd Parabolic Flight Campaign). 
In parallel to those experiments, we provide a quanti- 
tative explanation of the observed cavitation patterns 
through numerical simulations of dissipative shocks in- 
side spheres, and derive an analytic model to predict 
the location of shock-induced cavitation. Thereby we 
demonstrate that shock-induced cavitation in isolated 
volumes is a universal phenomenon. Finally, we discuss 
a potential implication of cavitation caused by trapped 
shocks for drop erosion of solid surfaces. 

Our microgravity experiment can be seen as an ideal 
laboratory for the study of shock dynamics inside sta- 
ble drops. The setup [details in Ref. Il5| can pro- 
duce a spherical drop of demineralized water (diam- 
eter D = 16 — 26 mm). This drop is smoothly expelled 
through an injector tube, which also serves as a perma- 
nent attach-point for the drop (see Fig.[2l left). A mov- 
able pair of electrodes penetrating the drop from the 
top in Fig. [2] releases a fast (10 ns) discharge at a spe- 
cific location within the drop. This discharge forms a 
sup ersonically expanding point-plasma [physics in Ref. 
[l6| , which generates a spherical shock wave and a va- 
por bubble pJTj, called the "primary cavitation bub- 














FIG. 1: Examples of cavitation (arrows) produced by re- 
flected shock waves, (a) Liquid drop (88% water, 12% gela- 
tine, diameter D — 10 mm) impacting on a hard solid sur- 
face at V — 110 ms~^ [reprint from Ref. J3]; (b) Water jet 
(D = 6 mm) shocked by a laser pulse (energy = 30 mJ) on 
the left side [setup explained in Ref. [3]; (c) Water jet {D = 
22 mm) impacted by 9 mm-projectile v — 200 ms~^ [col- 
laboration with the Swiss army] . Dashed lines indicate the 
position, where the cavitation should occur according to 
eq. (|2}, corrected for optical refraction 15]. 



ble" . We here focus on the shock wave, while regard- 
ing the primary cavitation bubble as a welcome side- 
effect to estimate the energy of the shock wave. In fact, 
early studies of laser-induced point-plasmas [15, Ts] ev- 
idenced that the shock energy approximately equals 
the bubble energy = An/Z- r^g^^ • (Poo —Pv)^ where 
^max is the maximal bubble radius, Poo = 80kPa is 
the static water pressure (= ambient pressure in the 
aircraft cabin), and = 3.2 kPa is the water vapour 
pressure at the working temperature T ^ 25 °C. Here 
we study cases in the range rmax = 2 — 4 mm, hence 
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~ ^ 3 — 20 mJ. The fast phenomena produced 
by the shock are recorded using a high-speed camera 
(Photron Ultima APX) at up to 120 000 frames/sec 
and 50/im spatial resolution. 

The data is acquired in 36 microgravity cycles, each 
of 20 s duration. Each cycle allows us to generate a 
single water drop and to release one shock wave within 
this drop. The free parameters defining such a cy- 
cle are the drop diameter the shock energy E^, 
and the location of the shock center. The latter is 
expressed by the eccentricity e = s/{D/2) G [0,1], 
where s is the distance between the drop center and 
the shock center. The results presented in this letter 
rely on the 18 parameter configurations corresponding 
to all possible combinations of G {16 mm, 22 mm}, 
Es e {3mJ,10mJ,20mJ}, and e G {0.2,0.4,0.6}. 
Each configuration is repeated twice to verify the re- 
peatability of the data. 

All 36 high-speed visualizations exhibit the same 
dominant feature: About 10 /is after the shock gen- 
eration thousands of submillimetric bubbles appear 
synchronously in the hemisphere opposite the origin 
of the shock (Fig. [2] left). The fastest visualizations 
(120 000 frames/sec) uncover that these micro-bubbles 
grow and collapse, their longest life-times lying around 
50 /Lis. This transient behavior and the coherent forma- 
tion of the bubbles at the instant of the shock transi- 
tion disclose that the bubbles are a form of shock wave- 
induced cavitation p]-H, fiol, [2Qj - a phenomenon known 
as "secondary cavitation" in other situations. The 
largest bubbles have diameters of d ^ 0.5mm, which in 
the Rayleigh-model \2^\ of spherical cavitation bubbles 
implies life-times of = 0.915 d [p/(poo — Pv)]^^'^ ^ 
50 /is, consistent with our observations. 

The discrepancy between the high bubble-density 
seen in Fig. [2] (left) and the faint traces detected by oth- 
ers ^ can be explained by the particular confinement 
of our shock. The free drop surface causes an elastic 
reflection of a pressure wave; hence the shock bounces 
back and forth while successively dissipating its energy 
to shock-induced cavities. To check if the full shock en- 
ergy is converted into cavitation bubbles, we estimate 
the total volume V of these bubbles from their radii 
(corrected for optical refraction by using our optical 
model in Ref . ITHf) . The implied energy E = V-{poo—Pv) 
is systematically consistent (within 20% measurement 
uncertainties) with the original shock energy E^. 

What is the role of gas and impurities contained 
in the water? Shock-induced cavitation bubbles arise 
when the shock pressure excites pre-existing nuclei [l| . 
A priori, one could therefore expect that the amount 
of shock-induced cavitation inside a shocked drop de- 
pends on the initial nuclei content. However, in the 
particular case where the shock keeps bouncing off the 
free surface until its entire energy is converted into bub- 
bles, the nuclei content cannot play a major role for the 
amount of shock-induced cavitation. Our observations 
of strong shock-induced cavitation in very clean water 
confirm this peculiar feature of confined shocks. 

Depending on the experimental parameters, the 
cloud of shock-induced bubbles appears in different in- 
tensities and geometries (Fig. [2]). A systematic inves- 
tigation of all samples uncovers that the intensity (size 
and number) of bubbles varies with the shock energy 
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FIG. 2: Spherical water drop (D = 22 mm) 80 /is after 
the generation of a spherical shock at the location labelled 
"cavity" (since here the primary cavitation bubble forms in 
the case of a plasma-driven shock). Left images show the 
observed shock-induced cavitation, while the right images 
are simulated counter-parts; the colors range from lowest 
(blue) to highest (red) energy density. The eccentricities of 
the shock origins are e = 0.4 (upper) and e = 0.6 (lower), 
though they appear to be larger due to optical refraction 
[refraction model in Ref. 15]. Dashed crosses indicate where 
the strongest cavitation should occur according to eq. ([2|). 



Es and drop diameter D. In return, the geometry of 
the bubble-cloud exclusively depends on the eccentric- 
ity e, as will be explained hereafter. 

To understand the patterns formed by shock-induced 
cavitation bubbles, we implement a new simulation 
scheme of dissipative shocks inside spherical drops 
based on a smoothed particle approach. In our model, 
the shock is considered as a spherical shell that isotrop- 
ically expands from a single point, specified by its ec- 
centricity e. This shell is represented by = 10^ par- 
ticles, which initially propagate on straight lines at the 
speed of sound c = 1500 ms~^. When reaching the 
drop surface they are reflected elastically. Each par- 
ticle carries a shock energy Ep{t) oc (^p)^, where 6p 
is the pressure fluctuation relative to the equilibrium 
state. (Note that Sp changes its sign when reflecting at 
the free surface, while the energy remains unchanged.) 
As these "shock quanta" propagate across the liquid 
they form cavitation via evaporation in the tensile case 
{5p < 0), or through forcing the collapse and subse- 
quent rebound of existing nuclei in the compressive 
case {dp > 0). In both cases, the generation of cavita- 
tion bubbles decreases the energy of the shock by de- 
creasing \Sp\. We assume that this decrease in energy 
happens at a constant fraction per unit time (indepen- 
dently of the sign of Sp) , as typical for most dissipative 
processes. Thus, Ep{t) = Eg/N - ex.-p{—t/r)^ where r is 
the time, in which the shock dissipates (1 — 1/e) ^ 63% 
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FIG. 3: Circular wave reflected inside a circle. The thick 
line highlights a selected path. The ensemble of all re- 
flected paths defines the dashed envelope, called "catacaus- 
tic" . This envelope has a peak point /i that depends on the 
wave origin e via eq. ([2|). Here, e = 0.4 as in Fig. [2] top. 



of its energy. Here we choose r = 20 /is since the 
shock-induced bubbles in Fig. [2] (left) reach an inte- 
grated potential energy, measured through the bubble 
sizes, of (0.6 ± 0.2)Es within 20 /is. The simulation 
progresses at a discrete time step dt = 10~^ s, cho- 
sen sufficiently small to obtain converging results. In 
each step and for each particle, the dissipated energy 
\dEp{t)\ = —Epdt is transcribed to the liquid at the lo- 
cation of the particle using a 3-d Gaussian smoothing 
kernel with a variance equal to the mean separation 
between neighboring particles. To tackle the energy 
acquired at different locations in the drop, the latter is 
discretized on a regular mesh of 256^ ^ 1.7 • 10^ cubic 
cells. 

The output of these numerical simulations is a 3-d 
density map of the energy dissipated by a spherical 
shock trapped inside a spherical drop. To compare 
these maps against the 2-d images of shock-induced 
cavitation inside drops, the simulated maps were pro- 
jected onto a plane, while accounting for the optical 
refraction at the water surface [according to our model 
in Ref. [l5[. Fig. [2] compares the simulation against 
the observations for the two eccentricities e = 0.4, 0.6. 
The patterns of the simulated energy density and the 
observed shock-induced cavitation exhibit remarkable 
similarities: in both cases (i) the highest densities 
appear opposite the origin of the shock wave at a 
comparable eccentricity, and (ii) for high eccentricities 
(e > 0.5) the bubble-patterns spread out in two wings. 
This match between simulation and experiment con- 
firms that the patterns of shock-induced cavitation can 
be understood in terms of reflected dissipative waves. 

Can we then analytically understand the observed 
and simulated patterns of shock-induced cavitation? 
Due to axial symmetry, the 3-d geometry of a spheri- 
cal wave inside a sphere can be reduced to a 2-d cir- 
cular wave reflected inside a circle (Fig. [3]). Unlike 
ellipses circles yield no focal point for an eccentric cir- 
cular wave. Instead, the reflected rays are concentrated 
within a zone enveloped by the "catacaustic" (dashed 
line in Fig. [3|), defined as the location where reflected 
rays intersect. On this catacaustic, we expect a high 
density of the reflected shock. This model can be ex- 
tended to 3-d by performing a rotation about the axis 
of symmetry (horizontal axis in Fig. [3|). Doing so, we 
introduce an additional 1/r-factor in the shock den- 
sity, where r is the distance from the axis of symme- 



try. Thus, the location of the highest density of shock- 
induced cavitation is defined by the intersection of the 
catacaustic and the axis of symmetry. This intersec- 
tion lies at a position fi (see Fig. [3]) that depends on 
the wave origin e G [0, 1]. To express /i as a function of 
e, consider a single ray emitted at an angle a G (0, 27r) 
and crossing the center line at the position x after its 
reflection (thick line in Fig. [3]). From the law-of-sines 
sin(/3/2) = e sina and sin(/3/2) = —x sin 7, and triv- 
ially a + /3 + 7 = TT, which solve to 



sin [a + 2 arcsin(e sin a)] 

Eq. ([1]) is meaningful if x G [—1,0], while otherwise 
the ray is reflected more than once before crossing 
the center line. For any e G [0,1], the maximum 
fi = max{x(e, a) G [—1,0]} is reached as a ^ 0, thus 

M = -e/(2e+l). (2) 

In conclusion, eq. (|2]) specifies the location of strongest 
shock-induced cavitation. This analytic prediction 
provides an excellent fit to the microgravity data as 
well as to the ground-based observations (see dashed 
crosses and lines in Figs. [T] and [2]). 

Motivated by the observational, numerical, and an- 
alytical evidence for shock-induced cavitation in drop 
and jets, we finally wonder about the potentially ero- 
sive implications of this cavitation. In fact, in Fig. [it 
the impact of the projectile onto the liquid jet is fast 
enough for the cavitation bubbles to survive until the 
projectile reaches them. Since cavitation can erode 
nearby surfaces [22[, this observation suggests that 
fast Hquid-solid collisions may cause erosion via shock- 
induced cavitation at the far-side of the impact point. 
To calculate the required impact velocity for a drop 
with diameter D (e. g. a rain drop), we assume that a 
spherical shock wave is emitted at the initial contact 
point. Hence, the strongest shock-induced cavitation 
is located at a distance 2D/3 from the contact point 
[eq. ([2]) with e = 1]. The reflected shock reaches this 
focus at a time tc ~ 41)/ (3c) after the impact, where c 
is the speed of sound, while the impacting solid itself 
reaches the same point at a time ~ 2D/{3v), where 
V is the impact velocity. For erosion to take place, we 
therefore require tc + ^R > ^v, where is the life-time 
of the bubbles [21|. Hence, a lower impact velocity 
limit for this type of erosion is 

Eq. (|3]) is valid for 2v < c, since otherwise the solid 
reaches the shock focus before the reflected shock it- 
self. The characteristic bubble life-time depends on 
the liquid and shock parameters, although our cases 
(Figs, m [T]3,c) all yield average life times of order ^r = 
20 /is (with the highest values reaching ^r = 50 /is). 
Adopting ^R = 20 /is for the case of a typical rain drop 
{D = 3 — 4 mm, c ~ 1500 ms~^) eq. ([3j) then gives a 
lower velocity limit of v ~ 100 ms~^, which is roughly 
an order of magnitude above the free fall velocity of 
rain. This reveals that the mechanism of cavitation 
erosion described in this section will only be active in 
particular cases, such as aircrafts and missiles eroded 



FIG. 4: Four mechanism of erosion by an impacting liquid drop (see main text). The last mechanism is proposed in this 
letter. The conditions for these mechanisms differ; hence only some of them may affect a particular impact. 



by rain j23|, and Pelton turbine blades eroded by 
high-speed (up to 200 ms~^) droplets and jets [25]. 

To complete the picture, Fig. |4] presents a synthetic 
view of the most important mechanisms of erosion that 
are known to occur during the impact of a liquid drop, 
including those addressed in the past (damages by such 
erosion are discussed in Refs. [231 [23, l26| ) : 

(a) Hammer pressure: On initial contact, erosion can 
result from the high compression-pressure of up to 
p = 3pcv for rigid solids [13], which exists while the 
contact edge expands supersonically. 

(b) Jetting: When the contact edge becomes subsonic, 
the pressure-shock detaches from the contact edge and 
high-speed jets emerge from the latter [l3|, causing 
shear erosion \27] (e. g. crater rims formed by hyper- 
sonic impacts of atoms [28] and meteorites [29]). 

(c) Near-side cavitation: A lateral pressure shock trav- 
elling from the contact edge to the axis of symmetry 



can produce cavitation next to the initial contact point 
[sol, lllj , which may cause point-like erosion [32] . 
(d) Far-side cavitation: New mechanism suggested 
here, which occurs via cavitation caused by reflected 
shocks opposite the impact point (see also simulations 
of low pressure at the shock- focus [33l-[35|). 

We emphasize that a particular drop impact may 
only involve some of these erosive mechanisms. A sys- 
tematic study of the applicability of each mechanism 
as a function of dynamical and geometrical parameters, 
as well as material properties promises an interesting 
road for forthcoming research. 
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